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We present a formalism for rigorous calculations of cross sections for inelastic and reactive colli- 
sions of ultracold atoms and molecules confined by laser fields in quasi-2D geometry. Our results 
show that the elastic-to-inelastic ratios of collision cross sections are enhanced in the presence of a 
laser confinement and that the threshold energy dependence of the collision cross sections can be 
tuned by varying the confinement strength and extern al mag netic fields. The enhancement of the 
elastic-to-inelastic ratios is inversely proportional to •^/e/HuJo, where e is the kinetic energy and ujq 
is the oscillation frequency of the trapped particles in the confinement potential. 



Atomic ensembles cooled to ultracold temperatures 
can be confined by optical forces of counterpropagating 
laser beams to form a periodic lattice structure. Optical 
lattices can be used to produce low-dimensional quantum 
gases by confining the motion of ultracold particles in 
one or two dimensions pQ. Ultracold atoms can be com- 
bined to form ultracold molecules. Molecules confined in 
low dimensions may undergo inelastic and chemically re- 
active collisions, which suggests the possibility to study 
chemistry in restricted geometries. Several previous stud- 
ies showed that collision dynamics of ultracold molecules 
restricted to move in two dimensions is different from 
scattering processes in an unconfined three-dimensional 
(3D) ultracold gas |3] H E] . The effect of the confining 
laser force on inelastic and reactive collisions of molecules 
in an optical lattice, however, remains unknown. Ultra- 
cold atoms and molecules in restricted geometries may 
be used for quantum simulations of fundamental many- 
body systems [6j[7] and the development of novel schemes 
for quantum computation [8] [9] [10]. An analysis of in- 
elastic scattering in a quasi-2D trapped gas is necessary 
to understand the feasibility of the quantum simulation 
proposals |6j [7J H] [10] . Ultracold atoms and molecules 
trapped in a quasi-2D geometry can also be used as con- 
trollable model systems of excitons and exciton polari- 
tons in microcavity semiconductors jlll IT2l ITS"] . Stud- 
ies of inelastic interactions in quasi-2D ultracold gases 
may thus find applications for new research in chemi- 
cal physics, quantum condensed-matter physics, quan- 
tum optics of semiconductors and quantum information 
science. 

Collision dynamics of atoms and molecules at ultracold 
temperatures is determined by Wigner's threshold laws 
[2] H3], which give the energy dependence of the scatter- 
ing cross sections in the limit of low energies. The thresh- 
old laws change with the dimensionality of the system [2]. 
However, the interaction of molecules confined by a har- 
monic laser force in one dimension cannot be described 
as a collision process in two dimensions j3j2]. Molecules 
confined in a one-dimensional optical lattice move freely 
in two dimensions and oscillate harmonically in the third 
dimension, which is usually referred to as a quasi-2D ge- 
ometry. At the same time, the interaction forces between 



the colliding molecules are much stronger than the laser 
confinement. The reactive complex of molecules is there- 
fore unconstrained and the reaction process occurs in 3D. 
An inelastic collision or chemical reaction releases a lot 
of energy and accelerates the collision products, which 
are therefore free to move in 3D as well. The effect of the 
confining laser force is thus only to restrict the motion of 
molecules before the collision. Petrov and Shlyapnikov 
developed a theory of elastic collisions between atoms in 
a quasi-2D gas based on the renormalization of the scat- 
tering wave function [3J. Here, we extend their work to 
develop the formalism for quantum calculations of prob- 
abilities for inelastic and chemically reactive collisions 
of molecules confined in quasi-2D geometry. Our work 
leads to an important conclusion that the ratio of cross 
sections for elastic and inelastic collisions is enhanced in 
the presence of laser confinement. 

In a quasi-2D system, the strength of the laser con- 
finement can be described by the oscillation length Iq = 
\J hj fJWQ, where ljq is the frequency of the harmonic po- 
tential and fi is the reduced mass for the collision com- 
plex. The oscillation length of the confining potential is 
usually much larger than the characteristic radius r e of 
intermolecular interaction potentials [3]. Therefore, at 
short interparticle separations r < r e , the interaction be- 
tween collision partners is not affected by the confining 
potential and the wave function in the region of r be- 
tween r e and the characteristic de Broglie wavelength of 
the particles A e is proportional to the 3D s-wave scatter- 
ing wave function [3]. The theory of Petrov and Shlyap- 
nikov j3] relates the quasi-2D scattering wave function in 
the region r e <C r <C A e to the 3D wave function by a 
proportionality coefficient r]ipo(0), where (po is the wave 
function for the lowest energy oscillation in the confin- 
ing potential. In the present work, we assume that the 
temperature of the confined gas is much smaller than Hujo 
and molecules prepared in state a are trapped in a quasi- 
2D geometry. We consider collision processes that induce 
transitions from state a to another state a' and assume 
that any transition a — > a' results in loss of confinement. 
The indices a and a' are used to describe the internal 
energy as well as the chemical identity of the colliding 
particles, i.e. they specify the collision channels. The 
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probability of inelastic or chemically reactive collisions 
in 3D scattering is described by the elements S a >^ a of 
the scattering S- matrix. At r > r e , different collision 
channels are uncoupled. We express the s-wave compo- 
nent of the wave function for the confined channel a at 
r e « r < A £ as a regular single-channel wave function 
in 3D multiplied by r)ipo(0) 
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where k a is the wave number of the collision complex, 
(j> a is the eigenfunction of the Hamiltonian at r = oo and 
S a ^ a is the S-matrix element for elastic scattering in 
3D. Repeating the procedure of Petrov and Shlyapnikov 
[3] and using Eq. (1), we obtain the following expression 
for r\ in terms of the S'-matrix element 
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Here, e is the collision energy and w(e/2hoj ) is a complex 
function given by w(e/2fudo) = ln(BhuJo/ire) + vk |5], In 
the limit of k a — ► 0, the matrix element S a ^ a is related 
to the scattering length a of ultracold particles as a — 
{Sa^u 1) /2ik a . 

In the region r e <§; r <C A e , the wave function (1) for 
channel a can be generally written as 
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where £? a and SS a are the amplitudes of the incoming 
and outgoing scattering waves, and v a is a normalization 
constant [15] ■ At r e <C r <C A s , £/ a can be written as 
xA a , where A a is the amplitude of the incoming scat- 
tering wave for s-wave collisions in 3D. Comparing the 
coefficient in front of the term e~ %kaT in Eq. (1) with that 
in Eq. (3) and using the conventional form A a = i\pK jk a 
[16] . we obtain the coefficient \ 
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The confinement thus modifies the amplitude of the in- 
coming scattering wave in the 3D collision region r e <C 
r <C A e . Since the asymptotic motion of the collision 
products after a reactive or inelastic process is uncon- 
strained, a combination of the spherical Hankel functions 
and 3D spherical harmonics should be used to describe 
the wave function in the outgoing collision channels 
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where V is the rotational angular momentum (partial 
wave) of the collision complex in state a' and m\ is its 
projection on the quantization axis. The 3D wave func- 
tion after a collision (Vw) is related to the 3D wave func- 
tion before the collision (V>in) by the S'-matrix operator 



Vw = Sipin.. Therefore, the amplitudes of the outgoing 
scattering waves £$ a 'l'm' are related to the amplitude of 
the incoming wave s4 a by the S-matrix elements 
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Since the colliding particles are initially prepared only 
in state a, the scattered part of the wave function in the 
inelastic channels is given by 

V\ s nd = ~EEE v~/r~ x S a . Vm ' l ^ a wX (7) 

a'^a V m'. 
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The scattering amplitudes for inelastic collisions are 
defined as follows 



Cel= V cc' 2 fa'^ a 



Comparing Eq. (7) with Eq. (8), we obtain 
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which yields the integral cross section for inelastic or re- 
active scattering in a quasi-2D gas 
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where ij is given by Eq. (2) and (po(0) = (-71)*- 

Eq. ( 10 1 shows that the cross sections for inelastic or 



reactive collisions in a quasi-2D gas depend on the 3D 
scattering length of the colliding particles in state a as 
well as the confinement strength. In order to illustrate 
the effect of these parameters on inelastic scattering, we 
present in Fig. 1 the results of rigorous calculations for 
collisions of 87 Rb atoms in the my = state with 6 Li 
atoms in the my = — | state, leading to Zeeman relax- 
ation in a magnetic field. The quantum numbers rrif 
denote the projections of the total angular momenta of 
the atoms on the magnetic field axis. The calculations are 
based on accurate interaction potentials for the 6 Li— 87 Rb 
molecule [T7]- The scattering length a of the 6 Li— 87 Rb 
system is tuned by varying an external magnetic field 
near the Feshbach resonance at 1104.9 G. According to 
Wigner's threshold laws, the cross sections for inelastic 
transitions in the limit of low collision energy e vary as 
- ^ in 3D [H] and as - ^l n 2 1 in 2D [5j. Fig. 1 
shows that the energy dependence of the cross sections 
for inelastic scattering in a quasi-2D gas resembles the 3D 
threshold law if \a\/lo << 1. It becomes similar to the 
energy dependence in 2D when \a\/lo > 1. This suggests 
that the threshold energy dependence of inelastic cross 
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FIG. 1: The threshold energy dependence of cross sections 
for inelastic relaxation in s-wave collisions of 6 Li with 87 Rb: 
filled circles— purely 2D geometry; filled squares— 3D scat- 
tering cross section reduced by a factor of 4 x 10 4 ; open 
circles— quasi-2D with \a\/lo > 1 cross section reduced by a 
factor of 30; open squares— quasi-2D with \a\/lo « 1. The 
initial states are mf Li = -| and mf Rb = 0. 



FIG. 2: Cross sections for s-wave inelastic collisions of Li and 
87 Rb atoms in 3D (solid curve) and quasi-2D scattering with a 
weak confinement (lo = 10 4 bohr - dotted curve) and a strong 
confinement (lo = 10 3 bohr - dot-dashed curve) as functions 
of the magnetic field. The inset shows the ratio of the cross 
sections for inelastic collisions in quasi-2D with lo = 10 4 and 
3D. The collision energy is 10 _8 cm _1 . The initial states are 
rrif Li = - ~ and mf Rb = 0. 



sections in quasi-2D systems can be tuned by varying 
the ratio |a|/7o. 

Fig. 1 also demonstrates that the laser confinement 
reduces the magnitude of the inelastic cross sections and 
that the suppression is more significant for the larger 
value of |a|/7 . In order to examine the dependence of 
the suppression on the scattering length, we present in 
Fig. 2 the cross sections for inelastic Zeeman relaxation 
in 3D and quasi-2D collisions of 6 Li and 87 Rb atoms as 
functions of the magnetic field varying through the reso- 
nance. The position of the Feshbach resonance is shifted 
by the confinement in agreement with previous calcula- 
tions of elastic cross sections [3] [18] [19]. The effect of the 
confinement is enhanced near the resonance. 

The inelastic cross section is suppressed because in 
the limit of zero collision energy it must smoothly ap- 
proach the threshold energy dependence for scattering 
in a purely 2D geometry [5]. In order to quantify 
the suppression of inelastic scattering, it is necessary 
to consider the ratio of cross sections for elastic and 
inelastic collisions. The elastic-to- inelastic ratio is of 
paramount importance for experiments with ultracold 
atoms and molecules. Elastic collisions determine the 
macroscopic dynamics of quantum gases. Inelastic and 
chemically reactive collisions destroy ultracold atoms and 
molecules. The ratio of cross sections for inelastic and 
elastic scattering in a quasi-2D gas can be written as 

CT quasi-2D^quasi-2D = ^D/^aD where ^SD denQte th(J 

cross sections in an unconfined 3D gas. Using Eqs. (2) 
and (18) of Ref. [3] and Eq. (10 1 derived here, we find 
that 7 = k a lo/\^K = \j2e/irhujQ. Because e/hujQ <C 1 in 
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FIG. 3: The ratios of inelastic and elastic cross sections for 
s-wave collisions of 6 Li and 87 Rb atoms as functions of lo for 
\a\ = 13.58 Bohr (B = 200G) (circles) and \a\ = 1704.43 Bohr 
(B — 1104. 9G) (squares). The collision energy is 10 _8 cm _1 . 
The initial states are mf Li = -i and mf Rb = 0. 



a quasi-2D gas, the ratio of elastic and inelastic cross sec- 
tions must always be enhanced under laser confinement 
and must increase as the laser confinement increases. The 
degree of the enhancement is given quantitatively by the 
equation above. This general result is illustrated by a 
numerical calculation in Fig. 3. 

The formalism presented in this paper can be applied 
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to describe chemical reactions in an ultracold molecular 
gas under laser confinement. The index a' in Eq. (10) 
must then include outgoing channels in different chemical 
reaction arrangements. To explore the effect of laser con- 
finement on chemical interactions of ultracold molecules, 
we consider an illustrative example of the reaction 7 Li 
+ 6 Li 2 (w = 0,N = 1) -^ 7 Li 6 Li + 6 Li. The cross sec- 
tions for elastic scattering and chemical rearrangement 
transitions in 7 Li + 6 Li 2 collisions have been calculated 
by Cvitas et al [20 . Using their results and Eq. (10) of 
this paper, we evaluate the cross sections for the chemi- 
cal reaction in a confined gas (Fig. 4). We note that the 
results presented in this paper apply to a gas of parti- 
cles trapped in the ground state of the laser confinement 
potential. The limit of 3D scattering therefore cannot 
be obtained from our results simply by increasing Iq. As 
described by Petrov and Shlyapnikov [3] , the limit of 3D 
scattering is obtained by heating the system so that the 
particles populate a manifold of states in the confinement 
potential. 

As demonstrated by Cvitas et al [20], rigorous quan- 
tum calculations of cross sections for ultracold reactive 
scattering in 3D are computationally very demanding, 
though not impossible. Breaking the symmetry of space 
by applying an external field increases the complexity of 
the scattering problem to a great extent and converged 
reactive scattering calculations in the presence of exter- 
nal fields are at present impossible (21]. The theory pre- 
sented here makes the analysis of reactive scattering of 
molecules in a quasi-2D gas feasible. 
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FIG. 4: The ratios of cross sections for elastic scattering (cir- 
cles) and chemical reaction (triangles) in quasi-2D and 3D 
as functions of the confinement strength for 7 Li + 6 Li2(v = 
0, N = 1) collisions. The collision energy is 10 _8 cm _1 . 



In summary, we have developed a formalism for rig- 
orous calculations of cross sections for inelastic and re- 
active collisions of ultracold atoms and molecules con- 



fined in quasi-2D geometry. The approach provides ex- 
pressions for the inelastic and reactive scattering cross 
sections in terms of the S-matrix elements for collisions 
in 3D and the laser confinement parameters. We use the 
method to elucidate the general features of inelastic scat- 
tering and chemical reactions in ultracold quasi-2D gases 
of atoms and molecules. We have found that the cross 
sections for inelastic and chemically reactive collisions are 
suppressed by the confinement forces. This suppression 
is generally more significant than the effect of the laser 
confinement on the probability of elastic scattering. The 
elastic-to-inelastic collision ratios are therefore enhanced 
in the presence of a laser confinement. Our results sug- 
gest that applying laser confinement in one dimension 
may stabilize ultracold systems. We have found that the 
threshold energy dependence of cross sections for both 
elastic and inelastic collisions in quasi-2D gases depends 
on the scattering length of the collision partners in the 
confined state and the confinement strength so it can be 
tuned by varying the confinement forces and an external 
magnetic field. This can be used to study fundamen- 
tal physics of threshold collisions and simulate dynamics 
of exotic quantum systems such as exciton polaritons in 
crystalline organic microcavities jllj . 
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